SOME NONLINEAR SECOND ORDER EQUATION MODELLING 

ROCKET MOTION 



DOROTA BORS AND ROBERT STANCZY 

Abstract. In this paper, we consider a nonlinear second order equation modelling 
rocket motion in the gravitational field obstructed by the drag force. The proofs of the 
main results are based on topological fixed point approach. 

1. Introduction 
In this paper we consider the following second order differential equation 

(1.1) z (t) = a (t) (z (t) + a) 2 exp (-^\ + P (t) 

with suitably chosen functions a, (3 to be defined below by (11.91) and some constants a, H 
to be determined later and suitable initial or boundary conditions derived from some 
model describing the motion of the rocket, cf. [4]. Specifically, it can be derived from the 
one-dimensional first oder differential equation 

(1.2) m (f) v (*) = {w (t) - v (t)) m (t) + f ext (t) . 

This equation can be used to model a rocket motion with decreasing mass m = m(t) and 
under external force f ext (t). The rocket should be equipped with the engine that emits 
mass mi = mi (t) with velocities w = w (t) , t G [to, ^i] • Therefore, the difference between 
the mass of the rocket m (t) and the mass of the gas emitted at time t is constant and 
equal m which is some positive value, i.e. m(t) = m x (t) + m . In m we include all 
components of the rocket infrastructure like the engine, empty tank, guidance equipment, 
payload and the like. The function f ext = f ext (t) denotes all external forces that stimulate 
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the rocket motion for example the gravity force or the drag force and in fact may depend 
also on the position and the velocity of the rocket, cf. (DQ) 

The expression mv(t) is the force acting on the rocket body which is given by the 
amount of mass expelled and the velocity of the mass relative to the rocket body, i.e. 
w(t) — v(t)rhi(t). The rocket moves forward since the rocket loses mass as it accelerates, 
so rhi(t) must be negative for any rocket. The term w — v^- is called the thrust of 
the rocket and can be interpreted as an additional force on the rocket due to the gas 
expulsion. The detailed derivation of the rocket equation one can find, among others, in 
[3], and 0. 

The equation (11.21) is so called the Mescerskii rocket equation, for details see [1], [3], where 
alos the case of more that just one emitted mass was considered. If we assume, the exter- 
nal force to be zero and a constant relative velocity of the emitted mass, i.e. f ext = and 
w (t)—v (t) =: c < for t G [<i, t 2 }, then the equation (11. 2ft has the form of the well-known 
Tsiolkovskii rocket equation. 

The corresponding solution to the Tsiolkovskii equation with the initial conditions m(t\) = 
m and v(t\) — v can be easily found after integration of (11.21) and reads as follows 

/x ~ , m (t) 
v (t) = v + c\n — — . 

m 

The Tsiolkovskii rocket equation in the gravitational field, i.e. f ext = —mg with the initial 
conditions m(ti) = m and v (ti) = v has a solution of the form 

, s / x , m (t) 

v(t)=v- git- h) + cln — —. 

m 

Assume that f ext = —mg — D with D = ^pv 2 ACn where Cn is a drag coefficient, A is 
the cross sectional area of the rocket, p is the air density that changes with altitude x, 
and can be approximated by p = po ex P ( — If) where H w 8000m is the so called "scale 
height" of the atmosphere, and finally po is the air density at sea level. In that case it 
is impossible to solve the equation explicitly. However, it is interesting to note that the 
effect of drag loss is usually quite small and is often reasonable to ignore it in the first 
approach. In practice, the drag force is approximately only about 2% of the gravity force. 
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Except for the rocket equation, the Tsiolkovskii equation can also be applied to the one 
photon decay of the excited nucleus. Moreover, it is possible to identify the Mescerskii 
equation with the bremsstrahlung equation where the electron is in the permanent decay. 
Many other decays can be investigated from the point of view of the Mescerskii equation, 
as one can see in [5]. 

In the paper we will consider the equation (jl.2p , with a given relative velocity of emitted 
masses q (t) = Vi (t) — v (t) for any i = 1, N and t G [ti,t 2 ] ■ Chemical rockets produce 
exhaust jets at velocities |q| from 2500 m/s to 4500 m/s. 

If we assume that the relative velocity is a given function c(t) — w (t) — v(t), we obtain 
the following equation 

m (t) v(t) = c (t) m (t) + f ext (t) . 

Let us put v (t) = x {t) where x it) denotes rocket position at time t. Furthermore, we 
assume that the external force depends not only on time t but also on the rocket position 
x (t) and its velocity v (t) according to gravitational and drag forces, i.e. 

f«rt/^ *U-U\ ..U\\ _^\- AC DP0„.2U\ ( X W 



r xt it) = f(t,x (t) , v (t)) = -m (t) g (t) exp 



H 



such that / : [t Q , tj x 1 x 1 f . 
Therefore, we get the equation 

(1.3) m{t)x{t) = c{t)m{t)-m{t)g-^^x 2 {t)expl X ^ 



2 v ' r \ H 
for any t G [to, h] . 

We provide our equation (jl.3p with the boundary conditions 

(1.4) x (t ) = x , x (ti) = xx, 

which guarantee that the rocket attains the given point in the required period of time. 
After a shift z (t) — x (t) — y (t) by a linear function 

(1.5) y(t) = a(t- t ) + x . 
where a = ^Eff an d equation (11. 3ft reads 

(1.6) z($ = c (t)±®_ g _ ACoPo (i (t) + a) 2 ^ ( z(t)+y(t) 



m(t) " 2 m(t) " V H 
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for any t G [toi^i] an d z satisfies the homogenous Dirichlet boundary conditions 

(1.7) z(t ) = z(t 1 ) = 0. 

Consequently, one can reduce considered equation to a aforementioned form 

(1.8) z{t) = a (t) (z (t) + af exp + P (t) 
where 

fi n\ U\ ac dPo ( y(t) 

(1.9) a(t) — -—exp' 



2m (t) "V H 

„ , , , s m (t) 

2. Integral formulation 

It should be noted that due to the presence of the derivative z the equation (11. 6p is not 
in a variational form. Therefore our approach differs from the one employed in [1] (for 
extension to elliptic problems see [2]) and makes use of integral formulation approach. 
Specifically the boundary value problem (II. 8ft with (j 1.7ft can be reformulated in the fol- 
lowing way 

ti 

(2.1) z(t) = jG(t,s)F(s,z(s),z(s))ds + b(t) 

to 

where the nonlinear term F is defined as 



F (s, z, z) — a (s) (z + a) 2 exp ( _ ^) 



with the abuse of notation, namely z denotes in the above formula the third independent 
variable and the function b is defined by 

h 

b(t) = Jg (t, s) /3 (s) ds. 

k 

Recall that a (t) and /3 (t) are defined by formulas ()1.9I) . The function G is the Green 
function, i.e. the continuous symmetric function satisfying the following conditions. 
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a. the homogeneous equation, i.e. for any t ^ s 

^(M) = 0, 

b. the boundary conditions, i.e. for any s G [to, £1] 

G(t ,s) = G(t u s) = 0, 

c. the jump condition, i.e. for any t G (t Q ,ti) 

,. dG , ^ dG , \ , 

hn l (*' s > ~ hm 7^7 (*' s ) = 1 • 

s^t+ at s-+t- at 

This function can be obtained from f ll.8f) by double integration and use of the the 
boundary conditions (II .TP and it reads 



(2.2) G(t, s ) 
with the derivative 



^{t- t ) {h-s),t< s, 
_i_( s _ to ) (t 1 -t),t>s. 



(2.3) fM= 5 f! tl "t t< " 

3. Main results 

We shall work in the space C 1 = C 1 ([to,ti] , M) of continuously different iable functions 
equipped with the natural norm 

||z|| = max {\z\ , |i|} 

where |z| = sup te r tot i \z (t) \ and \z\ = sup tg r t0jt i \z (t) \ for any z G C 1 . Next we can denote 
the right hand side of (12. ip as the value of the integral operator Sz, i.e. 

ti 

(3.1) Sz{t)= [G(t,s)F(s,z(s),z(s))ds + b{t). 



to 



Thus the existence of solution to (12. ip is reduced to finding a fixed point of the operator 
S, i.e. such azGC 1 that 

Sz = z. 
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First of all, note that the operator S is well-defined, continuous and compact. Indeed, for 
any z G C 1 the function Sz is continuous since the functions appearing under the integral 
sign ( 13. ip are uniformly continuous and (Sz)' is also continuous by the integrability of 
cf. ( I2.3P . The continuity of the operator S follows from the smoothness of F with 
respect both to z and z and boundness of the functions G and The compactness of 
the operator S requires the straightforward application of the Ascoli-Arzela theorem. 

3.1. Sign insensitive case. In this subsection we provide crude but immediate estimates 
yielding by the Schauder fixed point theorem the existence of a fixed point for the operator 
S. To this end we start with simple estimates of the following sup norms | • | 

(3.2) \Sz\ < G \ot\ {\z\ 2 + 2a \z\ + a) + \b\ 

2 



\(Sz)'\ < Gi \a\ (|i| + 2a \z\ + a) + \b\ 



where 

h 



G = sup / G (t, s) ds, 

te[t ,*i] J 



to 
h 



Gi = sup 

te[t ,ti] 



dG , , 



ds. 



to 



Using formula ( 12. 2 p one can easily derive the value of 

Go = — (ti — to) 2 



and 



G\ — (t\ — to) 2 • 



Theorem 3.1. The boundary value problem ( 11.7ft — ( 11.81) admits at least one solution 
provided (13.41) and (13.31) are satisfied. 

Proof. As announced we will make use of the Schauder fixed point theorem for the operator 
S. defined by ( 13.11) . We should prove that it maps some ball B (0, R) in C 1 into itself. Using 
the preliminary estimates (I3.2p and defining 

G 2 = \a\ max {G , Gi} = \a\ (ti — t Q ) 2 
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the invariance of the ball B (0, R) under the action of S can be reduced to the following 
condition 

G 2 (R 2 + 2aR + a) + \b\ < R 
which can be satisfied by some R > if 

A = (2aG 2 - l) 2 - AG 2 (G 2 a + |6|) > 0, 

i.e. 

( 33 ) | 6 | < 4a (a - 1) Gj - AaG 2 + 1 

4G 2 

and the value in the nominator is positive which is guaranteed provided 

□ 



Corollary 3.2. The conditions (13 Ah and (13.31) are satisfied if (3 is sufficiently small. 

Example 3.1. If A = 0, i.e. we neglect the drag force, accounting for negligable, in real 
world model, part, our assumptions are naturally satisfied. In such a context it is natural 
that one can attain the destination point in the desired period of time. In fact since the 
equation is linear one can derive the explicit formula for the solution. 

3.2. Sign sensitive case. In this section we provide more subtle a priori estimates tak- 
ing into regard the negative sign of some terms appearing in ( 11. 8ft and instead of using 
Schauder Theorem, we apply more subtle Schafer Theorem. First, we neglect negative 
term in (11.81) and integrate inequality 

z(t) < (3(t) 

using boundary conditions. Thus we get by ([1]) the bound for v and hence for z. Next we 
multiply by z both sides of (11.81) and integrate after dividing by (i + a) 2 ... to get better 
a priori bounds. 
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